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Abstract 

Superposition rules form a class of functions that describe general solutions of 
systems of first-order ordinary differential equations in terms of generic families of 
particular solutions and certain constants. In this work we extend this notion and 
other related ones to systems of higher-order differential equations and analyse their 
properties. Several results concerning the existence of various types of superposition 
rules for higher-order systems are proved and illustrated with examples extracted 
from the physics and mathematics literature. In particular, two new superposition 
rules for second- and third-order Kummer-Schwarz equations are derived. 
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1 Introduction 

The study of superposition rules can be traced back to the end of the 19th cen- 
tury, when Lie, Vessiot, and Guldberg [l]-[5] characterized and analysed the properties 
of systems of first-order differential equations admitting this property, the so-called Lie 
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systems [6]- [9]. Although the hnear superposition rule for homogeneous linear systems of 
first-order differential equations admits a natural analogue for homogeneous linear sys- 
tems of higher-order differential equations (HODEs), the generalisation of their nonlinear 
counterpart is not so evident and it has hardly been investigated so far |31 |T0] . 

Recently, the necessity of a theory of (linear and nonlinear) superposition rules for 
systems of HODEs became even more evident, as this concept repeatedly came up in the 
study of certain systems of second-order differential equations with multiple applications 
in physics and mathematics [T0]-|15]. 

In an attempt to fill in this gap of the mathematics literature, the present work 
aims to formalize the superposition rule notion for systems of HODEs and to analyse its 
properties. Since superposition rules for systems of second-order differential equations 
(SODEs) represent one of the most relevant types of superposition rules appearing in the 
literature, special attention is paid to this case. 

A notion of superposition rule for systems of SODEs was introduced in [10]. Neverthe- 
less, that work was more focused on the practical use of the concept than on studying its 
properties. That is why we start here by motivating this definition in detail and analysing 
some of its properties. 

The fundamental problem on the analysis of superposition rules for systems of HODEs 
is to find coordinate-free geometric conditions ensuring their existence. This problem, 
solved by the Lie-Scheffers theorem for systems of first-order differential equations, is 
here explicitly solved for systems of SODEs. Our new result provides not only a new 
insight into the study of superposition rules for SODEs, but also shows the existence 
of new and more powerful types of superposition rules for such equations. These new 
notions can be regarded as generalisations of other concepts already defined for systems 
of first-order differential equations (see ^). In addition, most of our achievements can 
be directly generalized to all systems of HODEs and they are also employed to review 
previous notions dedicated to the study of such systems, e.g. SODE Lie systems. 

Apart from their mathematical interest, our results are also relevant so as to study all 
physical systems and problems, like nonquadratic Hamiltonians or Berry phases (see [16] 
and references therein), related to differential equations admitting a superposition rule, 
such as second-order Riccati equations [TU] or Milne-Pinney equations [T5] . 

To highlight the interest of our methods, they are illustrated by the analysis of ex- 
amples extracted from the physics and mathematics literature. Special attention is paid 
to second- and third-order Kummer-Schwarz equations, whose mathematical interest is 
due, for instance, to their appearance in Rummer's problem, the study of Schwarzian 
derivatives, and other related topics [T71 \TE\ . Furthermore, Kummer-Schwarz equations 
occur in the analysis of non-stationary two body problems [T^l 120] and, via their relation 
to Riccati and Milne-Pinney equations [211 122], they can be employed to study several 
problems appearing in cosmology, quantum mechanics, and other branches of physics 
[T6l [T9| [2T| [22] . We here derive superposition rules for the analysis of such equations 
that provides us with several advantages with respect to previous methods of studying 
these, and other related, equations [211 [221 [23] • As a byproduct, we find a new property 
of Kummer-Schwarz equations: their dynamics is determined by a curve in a Lie algebra 
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of vector fields isomorphic to s[(2,M). 

The content of the paper is structured as follows. In Section 2 we describe some no- 
tions and results of the theory of Lie systems to be used throughout the paper. Section 3 
concerns the motivation and analysis of the definition of a superposition rule for SODEs 
as well as several particular types of it found in the literature. In Section 4 we provide a 
characterization of systems of SODEs admitting certain types of superposition rules and 
we describe a new kind of superposition rules for SODEs. In addition, several properties 
of superposition rules for SODEs are analysed. The relation of our new results and the 
so-called SODE Lie systems is studied in Section 5. The results of the previous sections 
lead to the definition and analysis, in Section 6, of a general notion of a superposition rule 
for systems of first- and higher-order differential equations. Subsequently, we illustrate in 
Sections 7 and 8 some of the theoretical results derived throughout our work by the inves- 
tigation of several remarkable HODEs. Finally, Section 9 summarizes our achievements 
and details some work to be accomplished in the future. 

2 Fundamentals on Lie systems 

We hereafter assume all geometrical objects and mappings, like vector fields or superpo- 
sition rules, to be real, smooth, and globally defined. In this way, we highlight the key 
points of our presentation by omitting the analysis of certain minor technical problems. 
For additional information, we refer to [91 [TO]. 

Definition 2.1. A superposition rule for a system of first-order ordinary differential equa- 
tions 

— =X^{t,x), ^ = l,...,n, (2.1) 
is a map $ : M™" x M"- -> M" of the form 

X = . . . , h,..., kn), (2.2) 

allowing us to write the general solution of system (12. ip as 

x{t) = <l>(a;(i)(t),...,a;(„)(t);A;i,...,fc„), (2.3) 

with (t), . . . , (t) being a 'generic' family of particular solutions and ki, . . . ,kn being 
the constants related to the initial conditions of each particular solution. 

Note 2.2. We shall not define rigorously what 'generic' means in the above definition, 
as it is not essential to our purposes and depends on the particular case. It shall be 
sufficient to bear in mind that, in the case of linear superposition rules for homogeneous 
linear systems of first-order differential equations, 'generic' means that the elements of 
the chosen finite family of particular solutions must be linearly independent. 

The uppermost achievement of the theory of Lie systems was obtained by Lie [1], 
who succeeded in characterizing systems of first-order differential equations that admit a 
superposition rule. 
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Theorem 2.3. (The Lie— Scheffers theorem) A system ^2.1\) admits a superposition 
rule Ii2. 2\) if and only if its right-hand side can he written as 



— = Z^it)e^ix) + ... + Zr{t)Cix), z = l,...,n, (2.4) 
so that the vector fields 

" d 

Xa{x) = J2Caix)g^. a = l,...,r, (2.5) 

i=l 

with r < m ■ n, span an r-dimensional real Lie algebra. 

The following definition and lemma, whose proof is a straightforward consequence 
of the Jacobi identity, notably simplify several statements and proofs of various results 
concerning the theory of Lie systems. 

Definition 2.4. Given a (finite or infinite) family A of vector fields on M", we denote by 
Lie(^) the smallest Lie algebra V of vector fields on M" containing A. 

Lemma 2.5. Given a family of vector fields A, the linear space Lie(,4) is spanned by the 
vector fields of 

A, [A, A], [A,[A,A]], [A, [A, [A, A]]],... 

where [A, B], with B = A, [A, A], . . . , denotes the set of Lie brackets between the elements 
of the families A and B of vector fields. 

Recall that if r : TR" — )■ R"- denotes the tangent bundle projection and 7T2 stands for 
the projection tt2 '■ {t,x) G R x R" 1— )■ a; G R", a time-dependent vector field X on R" is 
a map X : (t, x) G R x R" X{t,x) G TR" such that r o X = tts. Observe that every 
time-dependent vector field X on R" can be regarded as a family {Xt}t£R of vector fields 
on R'^, where Xf : a; G R" ^ Xt{x) = X(t, x) G T^R". 

Similarly to standard vector fields, time- dependent vector fields also admit integral 
curves 125]. We hereafter call an integral curve of X passing through [to, xq) G R x R" 
any integral curve : s G R )■ (t(s),7(s)) G R x R" of the one-dimensional distribution 
on R X R" spanned by the suspension of X, i.e. the vector field d/dt -\- X{t,x) [25] . 
satisfying that (to,Xo) G lm'y^°. 

From a modern geometric point of view, every system of first-order differential equa- 
tions of the form (12. ip is described by the unique time-dependent vector field on R", 
namely, X(t,x) = ^"=1 X*(t, x)(9/(9x% whose integral curves are (up to an appropriate 
reparametrisation) of the form {t,x{t)), with x{t) being a solution of system (12. ip . For 
simplicity, we use the symbol X to refer to both, a time-dependent vector field and the 
system of differential equations describing its integral curves. 

In such geometric terms, the Lie-Scheffers theorem states that a system X admits a 
superposition rule if and only if there exists a finite-dimensional Lie algebra of vector fields 
V, the so-called Vessiot-Guldberg Lie algebra, such that {XjjjgK C V . In consequence, 
X is a Lie system if and only if the Lie algebra 'Lie{{Xt}t&) is finite-dimensional. 
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The geometrical interpretation of superposition rules as well as one of the techniques 
for their determination is based on the notion of diagonal prolongation [9]. 

Definition 2.6. Given a time-dependent vector field X{t,x) = ^^^iX^{t^x)d/ dx^ on 
M", the time-dependent vector field X on of the form 



is called the diagonal prolongation to R"('"+^) of X . 

A method for determining superposition rules is briefly described as follows (see [HI ID] 
for details and examples). 

1. Take a basis Xi, . . . , X,, of a finite-dimensional Lie algebra fl2.5p associated with the 
Lie system under study. 

2. Choose the smallest positive integer m, so that the diagonal prolongations of the 
elements of the previous basis to (M")"^ are linearly independent at a generic point. 

3. Take global coordinates x^, . . . ,x" on W^. By defining this coordinate system on 
each copy of R"' within (M")'"^-'^, we get a coordinate system {x^f^^) M = 1, • • • , ^ = 
0, . . . , m} on (R'^)'""'"^. Obtain n functionally independent first-integrals Fi, . . . , F„ 
common to all diagonal prolongations Xi, . . . of Xi, . . . ,Xr to (R")'^"'"^ such 
that d{Fi, . . . , Fn)/d{x]^Qy . . . , x^g^) ^ 0. This can be performed, for instance, by 
means of the well-known method of characteristics. 

4. Assume the above first-integrals to take certain real constant values, i.e. Fi = ki for 
i = 1, . . . ,n. By means of these equations, calculate the expressions of the variables 
x^Q-) . . . , in terms of x^^-^, . . . , x"^-^, with a = 1, . . . , m, and fci, . . . , 

5. The obtained expressions give rise to a superposition rule in terms of any generic 
family of m particular solutions and the constants ki, . . . ,kn- 

Given two vector fields X and Y', we have that [X, y] = [X,y], i.e. the Lie bracket 
of two diagonal prolongations is a diagonal prolongation. Another, much less evident, 
property of diagonal prolongations is described in the following lemma, whose proof can 
be found in [9l Lemma 1]. 

Lemma 2.7. Consider a family of vector fields Xi, . . . ,Xr on R" whose diagonal prolon- 
gations to R"™ are linearly independent at a generic point. Then, given their diagonal pro- 



longations Xi,...,Xrto R"("^+i), a vector field X = ^^^^ fe^^a, with e (:7°°(R"("^+i)), 



is again a diagonal prolongation if and only if the functions are constant. 

It is worth noting that one can relate superposition rules to zero- curvature connections 



on a bundle pr : (x(o), . . . , X{rn)) G R"(™+^^ ^ . . . , X(rn)) e K"™ as follows (cf. [H]). 
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Proposition 2.8. Each superposition rule i[2.3\) for a system X is equivalent to a local 
n-codimensional foliation on R'^('"+^) whose leaves project, by pr, diffeomorphically onto 
M"™" and such that the vector fields {Xt}teR (^i"^ tangent to its leaves. 

The above result can be used to easily prove the following new result which is used 
posteriorly in order to analyse the existence of a particular class of superposition rules for 
systems of SODEs. 

Proposition 2.9. A family of Lie systems admits a common superposition rule if and 
only if they admit a common Vessiot-Guldberg Lie algebra. 

3 On the general definition of a superposition rule 
for SODEs 

To motivate the general definition of a superposition rule for SODEs, let us start by 
analysing a particular property of standard superposition rules. It is well known that 
every homogeneous linear system on M" of the form 

— = ^A^.(ty, . = l,...,n, (3.1) 



where the j are real t-dependent functions, admits its general solution x{t) to be written 
as 

X{t) = kiX(i){t) + . . . + knXl^n)(t), (3.2) 

with X(i)(t), . . . , X(^n)it) being a family of linearly independent particular solutions of (13. ip 
and /ci, . . . , /c„ a set of real constants. In other words, the system (13.11) admits a linear 
superposition rule. This leads to the existence of nonlinear systems admitting general 
superposition rules \8\. Indeed, every diffeomorphism : M" 9 a; t— ?■ 2; G transforms 
system (13. ip into 

= F\t,z), t = l,...,n, (3.3) 

at 

where the functions : M"'^^ — )■ M are generally nonlinear in the variables z^, . . . , z^, and, 
what is more important, whose general solution z{t) can be expressed (maybe nonlinearly) 

as 

zit) = 0(A;i0-^(2;(i)(t)) + . . . + M-^(2;(„)(t))), 

in terms of certain families of particular solutions Z(^i){t), . . . , Z(^n){t) of (13.31) and the con- 
stants ki, . . . ,kn. That is, since linearity depends on coordinate systems and the existence 
of superposition rules does not (recall the Lie-Scheffers theorem), the mere existence of lin- 
ear superposition rules for homogeneous linear systems of first-order differential equations 
leads to the existence of nonlinear systems admitting superposition rules. In addition, it 
is worth noting that not every system admitting a nonlinear superposition rule is of this 
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form. For instance, Riccati equations admit a superposition rule, but they cannot always 
be transformed diffeomorphically into linear homogeneous systems [26] . 

The aforementioned properties have an analogue for systems of SODEs. In fact, it can 
easily be proved that every homogeneous linear system of SODEs 

^ ^' - - ' • l,...,n, (3.4) 



dt 
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with A' j and j being any set of 2?t,^ time-dependent functions, admits its general 
solution to be written as 

X{t) = kiX(i){t) + . . . + k2nX(2n){t), (3.5) 

in terms of some arbitrary constants ki,...,k2n and a set of solutions {x(a)(t)|a = 
l,...,2n} such that the vectors {x[a)it),dx(a)(t)/dt) G TR" are linearly independent at 
every t G M. Now, a change of variables z = 0(x) transforms the above system into a 
(generally nonlinear) new one 

d'^z'' ^ ■ f dz\ , ^ 

W[t,z, — ], ^ = l,...,n, (3.6) 



rft2 \' ' dtj 

for certain functions : TR" x R — )■ R, and, moreover, such a change enables us, in view 
of (13. 5p . to write its general solution z(t) in the form 

Z{t) = <j){k^(j)-\z^^){t)) + ... + k2n<l)-\zi2n){t))), (3.7) 

in terms of a generic family of particular solutions 2;(i)(t), . . . ,Z(2n)(t) for (13. 6p and con- 
stants ki, . . . , k2n- Consequently, linear superposition rules for systems (13. 4p give rise to 
the existence of 'superposition rule-like' expressions for systems of SODEs. Expressions 
of this type frequently appear in the literature, e.g. in the study of linear inhomogeneous 
systems of SODEs. This suggests us the following definition that was proposed and briefly 
analysed in |10] and that includes the previous expressions as particular cases. 



Definition 3.1. A base-superposition rule for a system 

d'^x^ / dx 



dt- ^ V''"^'^'' ^ = 1'---'"' (3.^ 

is a map T : (M")"^ x R^" — )■ R" allowing us to write its general solution a;(t) as 

x{t) = T{x^i){t), . . . , x^rn){t); ki, . . . , k2n), (3.9) 

where X(i)(t), . . . ,X(^m)(t) is a generic family of particular solutions and ki, . . . , k2n are 
constants. 
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The above concept does not cover many other expressions found in the hterature for 
describing systems of SODEs [TOl |12[ [13], |15]. For instance, consider a Milne-Pinney 
equation 

^ = -a;^(t)x+^, (3.10) 

with X > and u{t) being any time-dependent real function [27j-[29]. This equation is 
relevant due to its applications in quantum mechanics, cosmology, Bose-Einstein conden- 
sates, and other physical topics [131 [El [16]. Recently, it was proved (see |13]) that its 
general solution can be written as 

1 /2 

x{t) = {hxl^it) +A;24)(t)±2 [Ai2[/34)(t)4)(t) -c(4)(t)+xj2)(t))]]'/'} , (3.11) 

by means of a generic pair X(i)(t), X(2)(t) of particular solutions, the function A12 = 
^12(^1, ^2, c. Is), the constant of motion 

/3=(^(iK.,(i)-^(i^„w)' + <^ 

and two constants ki,k2 related to initial conditions. Observe that expression (13. lip 
cannot be described by means of any base-superposition rule notion. Indeed, while ki, /c2 
take different values to describe the different particular solutions of (13.111) . the constant I3, 
whose value is fixed by the chosen particular solutions and their time-derivatives, does not 
appear in base-superposition rules. The same will happen for other new relevant expres- 
sions to be presented in this work. This motivates us to generalize the base-superposition 
rule as follows. 

Definition 3.2. A quasi-base superposition rule for a system of SODEs in M" of the form 
(13.81) is a function G : M™" x M'' x M^" — )■ R" allowing us to cast its general solution x{t) 
in the form 

x{t) = G(x(i)(t),...,X(„)(t),/i,...,/g;A;i,...,A;2n), (3.12) 

in terms of any generic family X(i)(t), . . . ,X(m)(t) of particular solutions of (13. 8p . a set 
of time-independent constants of motion Ii, . . . ,Iq, whose values are determined by the 
choice of the previous family and their derivatives with respect to the time, and a set of 
constants ki, . . . , k2n- 

Although almost every example of 'superposition-rule like' expression for SODEs is a 
particular instance of a quasi-base superposition rule, this notion still fails to cover sev- 
eral expressions found in the literature. That is the case of the very recently discovered 
expression for second-order Riccati equations, presented in [TU], which describes the gen- 
eral solution of such equations in terms of a generic family of particular solutions, their 
derivatives, and several constants. This motivates us to generalize the the concept of a 
quasi-base superposition rule as follows. 
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Definition 3.3. A system of second-order ordinary differential equations on R" given by 
(13. 8p admits a superposition rule if there exists a map T : (TM*^)"* x M^" — )• of the 
form 

X = T(X(i), f (1), . . . , Vi^rn)] h, ■ ■ ■ , hn) , (3.13) 

with {x(a),V(a)) £ Ta;^^jM" for a = 1, . . . , m, such that the general solution xit) of (13. Sp 
can be written as 

x{t) = T ^(t), . . . fci, • • • , h?j , (3.14) 

with X(i)(t), . . . , X(^rn)(t) being any generic family of m particular solutions of the system, 
and ki, . . . , /c2n being a set of constants related to the initial conditions of each particular 
solution. 

As every constant of motion involved in a quasi-base superposition rule can be con- 
sidered as a function on TM"*", quasi-base superposition rules can be easily regarded as 
a particular type of superposition rules for SODEs. Base-superposition rules can also 
be regarded as superposition rules that depend only on the base variables of TM"™' that 
justifies their name. 

Let us now turn to analysing several properties of superposition rules for systems of 
SODEs. Similarly to superposition rules for systems of first-order differential equations, 
expression fl3.14p cannot be applied for each family of m particular solutions. Recall 
that even in the simple case of a homogeneous linear system of SODEs, expression fl3.2p 
just remains valid for certain families of particular solutions. Consequently, in order to 
establish when a system (13.80 admits a superposition rule, it is essential to estabhsh what 
'generic' means in this new context. From now on, we say that (I3.14p is satisfied by a 
generic set of m particular solutions if there exists an open and dense subset U of (TR")™ 
such that expression (13.140 is valid for every family X(i)(t), . . . ,X(m)(t) that satisfies 

x,,m ^(0), • • • , x^o), %^(0)) e U c (TR")-. 

Every family of particular solutions satisfying the above condition is called a fundamental 
system of particular solutions of system (13. 8p . 

Superposition rules for systems of SODEs possess properties different from those of 
superposition rules for systems of first-order ones. Let us illustrate this fact by means 
of a particular remarkable difference. Consider again (12. ip as a Lie system admitting 
superposition rule (12. 3p . A time-reparametrisation r = T(t), with the inverse t = t{T), 
transforms this system into 

^ = ^F^(t(r),x), z = l,...,n. (3.15) 
dr dr 

As we assume (12. ip to be a Lie system, formula (12. 4p applies and the right-hand term of 
the above expression can be brought into the form 

uT^ (if 

^ = - (ZMrmiix) + ... + Zr{t{r))C{x)) , ^ = 1, . . . , n. (3.16) 
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Consequently, in view of the Lie-Scheffers theorem, the system fl3.15p becomes a Lie 
system. Moreover, as the general solution x{t) of (12 .ip and the general solution x(r) of 
fl3.15p satisfy x(t(r)) = x(r), then the superposition rule fl2.3p for (12. ip allows one to write 



XiT] 



$(x(i)(r),...,a;(^)(r);fci,...,fcn) 



in terms of a generic family a;(i) (r), . . . , (r) of particular solutions of system flS.lSp and 
ki, . . . ,kn- In summary. Lie's characterization of systems of first-order ordinary differential 
equations admitting a superposition rule is invariant under time-reparametrisations and 
Lie systems related in this way share a common superposition rule. Indeed, note that this 
follows trivially from the form of f l3.16p and Proposition 12.91 

The above property is no longer valid for superposition rules of systems of SODEs. 
Given a system of SODEs admitting a superposition rule, the systems obtained from it 
by time-reparametrisations do not necessarily possess the same superposition rule. For 
instance, consider a system of SODEs fl3.8p admitting a superposition rule (13.141) . A 
time-reparametrisation r = T(t), with the inverse t = ^(t), transforms (13. 8p into 

(fx^ (Pt dx^ dr f dt\^ ■ f , , dx dT\ 

+ ^M^M'^'^T^rh ^ = l,...,n, (3.17) 



dr'^ dr'^ dr dt \dT J \ ^ ^ dr dt 

whose general solution x(r) satisfies x{t) = x(t(r)), where x{t) is the general solution of 
(13. 8p . Hence, from the superposition rule (I3.14p . we get that x{t) can be expressed as 

<r) = T (a:(,)(r), • • • ,X(„)(r), ^(t(r))^(r); fc„ . . . , 

The above expression is not necessarily a superposition rule, as it may admit an explicit 
dependence on the new time variable r. A simple example illustrating this fact can be 
found in Section [71 

Obviously, we could have also required the superposition rule concept for systems 
of SODEs to be invariant under time-reparametrisations, but this would exclude several 
important examples like second-order Riccati or Milne-Pinney equations jlO[ [T5] . 



4 On the existence of superposition rules for SODEs 

The following theorem characterizes systems of SODEs admitting a superposition rule. 
We hereafter use canonical global coordinates (x^, . . . , x", v^, . . . , u") on TM". By defining 
this coordinate system on each copy of TM" within (TM")"^, we obtain a coordinate system 
{xl^ , I ^ = 1, . . . , n, a = 1, . . . , m} on (TR")™. 

Theorem 4.1. A mapping T : (TM")™ x M^" — > M" is a superposition rule for a system 
of SODEs (EJj if and only zf 
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1. the functions Uk : (TM")™ 3p^ Uk{p) = T(p; A;) G M", with k e 



p2n 



, are common 

solutions for the t -parametrized family of systems of PDEs on (TM") '" given by 

{X^^^)H + {Xt\ul = r{t, Uk, {X^^\uk), i = l,...,n, (4.1) 

where Uk = [u]., . . . ,u^) G M" and X^\x^l'^^ are the diagonal prolongations to 
(TR")'" of the time- dependent vector fields 



n / Q Q \ n 

i=l ^ ' i=l 



tF\t,x,v)^^ (4.2) 



and 

2. the map ^ : (TM")"^ x R^" TR" of the form 

^{p- k) = (ukip), [iX^j^\uk]ip)) G T„,(rtR" (4.3) 

gives rise to a family of bijections ipp : k E R^" ip{p\ k) G TR", with p being a 
generic point of (TR")™. 

Proof. Assume that the SODE Lie system (13. 8p has a superposition rule (12.31) . One can 
define the function : p G (TR")™- T{p; k) G R"", for each k G R^", which leads, 
for every fundamental system of solutions X(i){t), . . . , X(m)it) of (13. Sp . to a new particular 
solution of this system, 



xit) = Uk ( X(i)(t), ■ ■ ■ ,a:M(t), %^(t) ) . 



dt 



On the other hand. 



(4.4) 



dx 



dul 



du\ 



a=l i=l 



dx-'i 



(a) 



1 n, 



(4.5) 



where, for shortening the notation, we have denoted F^^^ = F^(t,X(^a),V(a)) and 

pW=(a^(i)W,^W,...,XM(t),^(t)). 
From the expression of X/) given in (14. 2p . it follows 



dx 



^(t) = [ix^-^)tuiMt)), I 



1 n. 



By differentiating expression (14. 5 p with respect to the time, we obtain 



(4.6) 



(4.7) 



it) 



a=l i=l 



EE 

, -^K^a) du 



d'^ul 



(a) 



+ 



dt dv^, 



+EE("i. 

1=1 i,«=i 



9F; 



(a) 



dx\ 



OF: 



(a) 



dul 



(a) 



(a) 
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If we compare the above expression with 





^^^^^fMTT-i — r lv,.t< (uy—, — — Vr,.l:<, 




we obtain 



a=l j=l ^^{a) a=l j,l=l \ ^^{a) ^^(a) / ^"(a) 

^'"^-(t) = iiX^^^)H + (Xf )),<)(p(t)), z = 1, . . . ,n. (4.8) 



As x{t) is a solution of system f l3.8p . and in view of expressions (14 .7^ and (14. Sp . it turns 
out that 

((^lr^)M + {Xt\^l){p{t)) = {t,uMt)), {{X^D^uMt))) . (4.9) 

Now, equation (14. 9 p holds for every fundamental system. This implies that, for each 
t G M, the above equation remains valid for a generic open and dense subset of (TM")'". 
Hence, 

iX^D^)tK + ixt^)tul = {t, n., iX^^\u,) , . = 1, . . . , n, 

for every t G M. Additionally, as the above procedure is still valid for every k G M^", every 
superposition rule provides us with a family of 2n-parametrized solutions Uk{-) = T{-;k) 
of the t-parametrized family of systems of PDEs (14.11) . 

Consider now a fundamental system X(i)(t), . . . ,X(m)(i) and denote p = p{0). For an 
arbitrary {xo,vo) G T^^IR", the theorem of the existence and uniqueness of solutions for 
systems of first-order differential equations shows that there exists a solution x{t) of system 
(13.81) with initial conditions a;(0) = Xq and dx/dt{0) = vq. In view of the properties of 
superposition rules, there exists a single k G M^"' such that x(t) = T (p(t); k) = Uk {pit))- 
Consequently, in view of expression (14. 7p . one has 



^0 = uUp), 

= [iX^D\<]iP), 



, _ dulipit)) 



-'0 



dt 



,n. 



t=o 

In other words, for a generic p G (TM")™, there exists a single k G M^" such that ip{p, k) = 
[xq, Vq). It follows that is a bijection that concludes the "if" part of our demonstration. 

Let us now prove that a map T : (TR'^)"* x M^*^ — )■ M" satisfying conditions (i) and 
(m) is a superposition rule for (13.81) . Consider any solution x(t) of (13.81) . Given a generic 
family of m particular solutions X(i)(t), . . . , a;(m)(t) of (13. 8p . condition (ii) ensures that 
there exists a unique k G M^" such that ip{p{0),k) = {x{0),dx/dt{0)), where we took 
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again p{t) to be of the form (14.61) . In view of condition (z), the function Uk{-) = T(-; k) 
is a solution for the family of systems of PDEs (14. ip . Defining now x{t) = Uk{p{t)) and 
using that expressions (14. 7p and (14. 8 p are valid again, we get 

^(t) = (t, x{t), ^{t)^ , ^ = 1, . . . , (4.10) 

That is, x{t) is a solution to (13.80 . Moreover, in view of condition (ii) and formula (14.70 . 
x(0) = x(0) and dx/dt{0) = dx/dt{0). Consequently, x(t) and x(t) are both solutions 
of (13. 8 p with the same initial conditions and they hence coincide. In summary, for every 
solution x(t) of system (13. 8 p and a generic family of m particular solutions, there exists a 
unique k G M^" such that a;(t) = T{p{t); k), so T is a superposition rule. □ 

Roughly speaking. Theorem 14.11 states that the existence of a superposition rule for a 
system of SODEs (13.81) is determined by the existence of an 'appropriate' 2n-parametric 
family of particular solutions of the family of systems of PDEs (14.11) . The interest of 
this result is obvious: it characterizes not only the existence of superposition rules for 
systems of SODEs, but also provides us with a tool, namely the family of systems (14. ip . 
to determine them. 

Note 4.2. Note that Xd and Xl are properly defined t-dependent vector fields over 
(TR"')™' and they maintain the form (14.20 for every coordinate system on (TM")"^ induced 
by a coordinate system on MJ^. 

Note 4.3. Denote by Sij a permutation of variables <H- xq), with i, j = 1, . . . ,m. As 
(14. ip and (14. 3 p are invariant under such permutations, it can be easily inferred that, if T 
is a superposition rule for (13.80 . then SijT is also, which provides an analogue for systems 
of SODEs of a known result about standard superposition rules [S]. 

Apart from the main result of Theorem 14.11 a careful analysis of its proof suggests 
us new types of superposition rules for systems of SODEs generalising previous notions 
used in the study of first-order differential equations [HI |30]. Indeed, given a particu- 
lar solution of the family of systems of PDEs (14. ip and a family of particular solutions 
X(i)(t), . . . ,X(^rn)(t) of system (13.80 . we can define 

The above expression has the same form as (14.40 . Following the calculations carried out 
in the "if" part of Theorem 14. H we obtain that the first and the second derivative of the 
above curve satisfy relations (14. 7p and (14. 8p . From here, as is a solution of (14. ip . it 
follows that x{t) is a new solution of (13. 8p . In other words, a particular solution of the 
systems of PDEs (14. ip allows us to generate new solutions of system (13. 8 p from any set 
of m particular solutions for this same system. This fact enables us to define a new type 
of superposition rule for systems of SODEs as follows. 
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Definition 4.4. A partial superposition rule for a system of SODEs fl3.8p is a mapping 
V : (TM")"' X RP ^ M", witli p < 2n, sucli tliat 

• For a generic set X(i)(t), . . . ,X(m){t) of particular solutions of system (13.81) . 



x{t) = V ( • • • ,a;M(t), ^^W; /ci, . . . , 



is a new solution of (13. Sp . 
• For a generic p G (TM")'", the map Vp :k & W ^ 'P{p\ k) G M" is an immersion. 

Obviously, for every fixed k = (fci, . . . , kp), the map Uj.{-) = V{ - ;k) is a solution of the 
system (14. ip . In view of this, it is easy to generalize Theorem 14.11 in order to characterize 
systems of SODEs admitting partial superposition rules. Moreover, the above notion 
extends to systems of SODEs the notion of partial superposition rule for systems of first- 
order differential equations defined in [S]. 

Let us now illustrate how the above results and definitions work. Consider the SODE 

and look for a superposition rule depending on a single particular solution. Following the 
terminology used in Theorem 14. H we have m = 1 (one particular solution) and n = 1 
(system defined on M). Consequently, the corresponding family of systems of PDEs (14. ip 
reads 

1^ + 2^(1)^ 7^ TT- + ^ irr + + 2^7^ = t\ (4.12) 



V 



'(1) <JX(^i)UV(^i) c-f(^) cx-(i) c-t;(i) 

whose common solutions, which do not depend on t, are solutions of the system 

d'^u ^ ^^0 2 1 

The solutions of the above system are of the form u = + ki, with ki being an 
arbitrary constant. Obviously, the family of systems (14.120 does not give rise to a two 
parametric family of solutions and (14. lip does not admit any superposition rule in terms 
of one particular solution. Nevertheless, it is interesting to point out that the solutions 
u = X(i) + ki exemplify that, for every particular solution X(i)(t) of (14.120 . the new 
function u{x^i)(t)) = X(i)(t) + A;i is a new solution of the system that gives rise to a partial 
superposition rule V : (x(i),f(i); ki) G TR x M t-)- + fci G M for equation (14.110 . 

Let us now turn to determining all possible superposition rules for (14.111) involving 
two particular solutions. So, we have m = 2, n = 1, and the family (14. ip reads 



( ^"^^ 2 ^ (t^ 2t ^ t 



2 
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Proceeding as before, we obtain that solutions of this t-parametrized family of PDEs are 
solutions of the system 

2 r\ 2 2 (-.o 

Ecu sr-^ O U sr^ O u 



2 r\'} 2 r^ 

Plugging the first equation of the above system into the others, we obtain that the above 
system is equivalent to 

du du du ^ du . .n "^^u 

whose solutions take the form u = — X(2))/i + + /2, with fi and /2 being 

two arbitrary functions depending on — V(2)- Now, by choosing appropriate one- 
parametric families of solutions of the above form, we can get partial superposition rules. 
For example, setting fi = k and /2 = 0, with A; G M, we obtain the family of solutions 
Uk = — X{2))k + which results in the partial superposition rule V{x{i),X{2y, k) = 
— X(2))k + which generates new particular solutions out of two known ones and 
one constant. Moreover, Theorem 14.11 shows that the determination of a superposition 
rule for system ( 14. lip amounts us to obtaining a two-parametric family of solutions u^kiM) 
of the above form such that condition {ii) of Theorem 14.11 holds. This can be done in 
several ways. For instance, by setting /i = ki and /2 = A;2, with fci, /c2 G IR, we obtain the 
superposition rule 

T(xi,fi,X2,f2; ki, k2) = M(fci,fc2)(a;(i),a;(2)) = ki{x(^i) - X(2)) + + /c2, 

and if we choose /i = — f(2)) and /2 = k2{v(^i) — f(2)), we arrive to 

T(xi, wi, X2, V2; ki, ^2) = M(fei,fc2)(a;(i), a;(2)) = (i) -^^(2))(a;(i) -a;(2)) +^2(^^(1) -^^(2))- 

Using our methods, we can easily derive the results of Table [H Special attention 
must be paid to the second example, illustrating that partial superposition rules may 
exist when superposition rules depending on the same number of particular solutions do 
not. In addition, this particular example is not a SODE Lie system (see definition in the 
next section), which was almost the only tool to study superposition rules for systems of 
SODEs so far. 

Theorem 14.11 characterizes systems of SODEs possessing a base-superposition rule as 
follows. 

Corollary 4.5. A mapping T : (M")"* x R^" — )■ is a base-superposition rule for a 
system of SODEs Ii3.8\) if and only if the following two conditions hold: 
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SODE 


Superposition rule 


Partial superposition rule 


IF ~ ^ 


+ ^2) 


+ ki 


dPx J. ( dx\^ 
df^ ~ ^\dt) 


nonexistent 


+ ki 


d''x J.2 dx 
dt^ ~ ^ dt 




kiX(i) 



Table 1: Superposition and partial superposition rules depending on a particular solution. 

1. The functions Uk '■ p & K"™" H- Uk{p) = T{p; k) G M", with k G M^", are solutions of 
the t-parametrized family of systems of PDEs on (TM")™' given by 



ix'r'')H + {Xr)tul = F\t, n, (Xr)t«.), ^ = i, . . . , n, 

where X^\x^2^'' are the diagonal prolongations to (TM")*" of the time- dependent 
vector fields 

" Q " 



(m)x 



2. The map (p : 



dx^ ^ dx^ 

1=1 4 = 1 

if{p; k) = {uk{p), [{x["'^)oUk]{p)) G 



gives rise to a family of bijections ipp : k E M?"' 1— )■ ip{j)\ k) G TR", with p being a 
generic point of (TM")™. 



5 Superposition rules and SODE Lie systems 

Recently, the theory of Lie systems was employed to obtain a few results about super- 
position rules for systems of SODEs [II]-[l5]. All these achievements were based on the 
notion of a SODE Lie system. We now describe this concept and provide several new 
results about the use of Lie systems to analyse different types of superposition rules for 
systems of SODEs. 



Definition 5.1. A system of second-order ordinary differential equations (13.81) is a SODE 
Lie system if the first-order system 



dx'' 

^ = F%x,v) 



i = l,...,n, 



(5.1) 



obtained from (13. 8p by adding the variables f * = dx^/dt, i = 1, . . . , n, is a Lie system. 
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The Lie-Scheffers theorem is an effective tool to determine whether a system fl3.8p is 
a SODE Lie system or not. Nevertheless, this method is based on analysing properties 
of the time-dependent vector field associated with the corresponding system (15. ip and it 
does not provide any straightforward information about the superposition rules for these 
systems. In order to overcome this drawback, we provide the following characterization 
of SODE Lie systems in terms of properties of superposition rules. 

Proposition 5.2. A system of SODEs liS. 8\) is a SODE Lie system if and only if it admits 
a superposition rule T : (TM")™ x M^" — )■ such that x|™'*T = 0, where Xl is given by 



Proof. Consider a second-order system of the form (13. Sp admitting a superposition rule 
f lS.lSp . The general solution x{t) of this system can be put in the form fl3.14p . Differen- 
tiating this expression with respect to the time, we obtain 



t w ^ E E wf- - E E ^ w (.2) 

a=l j=l (a) a=l j = l ^ ^ (a) 

where i = 1, . . . , n and p[t) is given by (14. 6p . Therefore, by defining 

f*(t,X(i),t;(i),...,X(^),t;(„)) = J2J2''U^-^ + J2J2^' {t,X(^a),Via)) (5.3) 

a=l j=l "^(a) a=l j=l "^{a) 

expressions fl3.14p and (15. 2p can be brought into the form 

X\t) = T • • • , X^m)it), ^(t); fci, . . . , fc^n 

^{t) = T (t, xa)(t), ^(t), . . . , x^t), ^(t); fc„ . . . , 

(5.4) 

Taking into account that the general solution {x{t),v{t)) of the first-order system (15. ip 
is obtained by adding the variables = dx^dt to system (13. 8p . we see that expressions 
([53D define a map $ : (t,p; fc) G M x (TR")™ x M^n ^ (x(p; fc), T(t,p; fc)) G TM" which 
allows us to write the general solution of this first-order system in terms of a generic set 
of particular solutions {x(a){t),V(a)(t)), with a = 1, . . . ,m. In view of expression (15. 3p . 

-^(t,X(i),t;(i),...,X(^),t;(„)) = 2^2^-g^{t,X(^a),Via)) —— = Xl . 

a=l j=l (a) 

Therefore, if x|™''T = 0, the mapping T and, in consequence, $ are time-independent. 
This shows that the function $ is a superposition rule for system (15. ip . which is therefore 
a Lie system. Hence, system (13. Sp is a SODE Lie system. 
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Let us assume now that system fl3.8p is a SODE Lie system, i.e. the first-order system 
(15. ip is a Lie system and there exists a superposition rule $ : (p; k) G TM"™ x M^" t— )■ 
(T(j9; k), $t,(p; fc)) G TM" such that its general solution {x{t),v{t)) can be written as 

X{t) = T(x(i)(t),t;(i)(t), . . . ,X(m)(t),f(m)(t); /^l, • • • , /^2n), 
V{t) = ^^{X(^i){t),V(i){t), . . . ,X(rn)it),V(^.ra){t); ki, . . . , Ajgn), 

where (x(a)(t), f(a)(t)), with a = 1, . . . , m, is a generic family of particular solutions of 
system Since dx(a){t)/dt = V(^a)it), the function T : TM"'" x M^n _^ enables us 

to write the general solution x{t) of system (13. 8p in the form 

X{t) = T (^X(i)(t), ^^(^), • • • , X^rn){t), —^(f)'^ h,---, k2n 

i.e. in terms of a generic family of its particular solutions X(i)(t), . . . , X(m,)(^), their deriva- 
tives, and a set of 2n-constants. In other words, our system of SODEs admits a super- 
position rule. Consequently, differentiating the above expression with respect to t, we 
obtain, in virtue of (15.31) . that $t, {pit)] k) = T (t,p{t); k) for a generic pit), which is given 
by (14. 6 p and constructed from a family of particular solutions X(i)(t), . . . ,X(m)(i)- Hence, 
, f (1) , . . . , X{rn) , f (m) ]k) =T (t, , , . . . , , f („) ; fc) and 

• • . fc) = (t, . . . , a;(m), k) = = 0, 

for i = 1, . . . , n. □ 

The above proposition improves the results of [lO], where it is only stated that SODE 
Lie systems admit superposition rules. Indeed, our new result also supplies additional 
information about such superposition rules, namely, x'ff^T = 0. This property is going 
to be used next to retrieve easily some previous results found in [TO] and various new 
ones. 

Proposition 5.3. Every system of SODEs Ii3.8\) admitting a quasi-base- superposition rule 
is a SODE Lie system. 

Proof. Assume that system (13. 8p admits a quasi-base-superposition rule. Let us now 
prove that this quasi-base superposition rule gives rise to a superposition rule for (13. 8 p 
such that X^^T = 0, which, in view of Proposition 15. 2[ proves that system (13. 8p is a 
SODE Lie system. 

The general solution x{t) of (13. 8p can be cast in the form (13.120 . As the functions 
Ij : TM"*" — )■ M, with j = 1, . . . ,g, are constant along the m-tuples {x(^a)it),dx(^a)it)/dt) 
obtained from m particular solutions X(a)(t) of system (13. 8p . i.e. 

Ui){t),^{t), . . . ,x^m){t),^^{t)\ = const., (5.5) 
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we obtain 



i=l a=l 



(a) \ / (a) 



where p{t) is given by expression (14.61) . The above holds for every generic family of 
particular solutions. Then, X^^Ij = 0, for j = Substituting expression (15. 5p 

into (I3.12p . it turns out that there exists a superposition rule T : (TM")'" x M.'^^ — )■ M" 
for system (13. 8p of the form T (p; k) = G{x(i), . . . , X(^m), h{p), ■ ■ ■ , Iq{p)i k), where p = 
(x(i), f (1), . . . , X(rn),V(m))- Indeed, in view of the definition of T and the properties of the 
quasi-base-superposition rule G, 

x{t) = G(x(i)(t), . . . , x^m){t), iMt)), . . . , iMt)); k) = r {p{ty, k) , 

where X(i)(t), . . . , X(m)(t) is any generic family of particular solutions of system (13.81) . 
Then, 

j=l i=l a=l ^ / [a) 3 j=i 3 

□ 

Corollary 5.4. Every system of SODEs admitting a base-superposition rule is a SODE 
Lie system. 

The implication of the above corollary cannot be reversed, i.e. not every SODE Lie 
system admits a base-superposition rule. Indeed, the following results can be easily used 
to prove the existence of SODE Lie systems admitting no base-superposition rule. 

Lemma 5.5. Given a system of SODEs h3. <^)) admitting a base-superposition rule, the 
systems 

d^x^ dr d'^t dx^ f dt\^ f , . dr dx\ , , 

+ {-r) F {ti^)^^^-ji-r] ^ ^ = l,...,n, (5.6) 



dr^ dt dr"^ dr \d.T J \ ^ ^ dt dr J 

with t = t(r) being any time-reparametrisation, are SODE Lie systems admitting a com- 
mon base-superposition rule. 

Proof. A time-reparametrisation t = t{T) maps system (13. 8p to a system of SODEs of 
the form with the general solution x{t) = x(t(r)). If T : M"™ x M^n ^ 

is a base-superposition rule for (13. 8p . then x(r) = T(x(i)(r), . . . , X(m)(T); /ci, . . . , A;2n), 
where X(i)(r), . . . , a:(m)(T) is any generic family of particular solutions of (15. 6p . Conse- 
quently, all the second-order differential equations of the family (15.60 admit a common 
base-superposition rule and, according to Corollary 15. 4[ all systems (15. 6p are SODE Lie 
systems. □ 
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Although only few SODE Lie systems admit them, base-superposition rules are the 
main superposition rules treated in the literature. The next proposition shows that SODE 
Lie systems must satisfy various restrictive conditions to admit a base-superposition rule. 

Proposition 5.6. Given a system of SODEs liS. 8\) admitting a base-superposition rule, 
the associated first- order system Ii5. 1\) is a Lie system related to a Vessiot-Guldberg Lie 
algebra containing the Liouville vector field of the tangent bundle TW^ and the vector 
fields 

^) = E rfA^^^^'^^^^i' ^' a^' P = 1. 2, . . . , (5.7) 

1=1 

where ti eR and \ eR'^ = R - {0}. 

Proof. In view of Lemma [5.51 time-reparametrisations r = T(t), with inverses t = t{T), 
transform system (13.81) into the family of SODE Lie systems (15.61) . whose associated first- 
order systems 

dr ' 

dv^ dr dH • f dt\^ ■ ( , . dr 



dr dt dr'^ \dT J \ ^ ^ dt 

with i = 1, . . . ,n, admit a common base-superposition rule. According to Proposition [231 
this implies that there exists a finite-dimensional Lie algebra of vector fields V containing 
all the vector fields 



with t = t(r) being any time-reparametrisation and r G R. In particular, if we take 
t{r) = r, we obtain that 



i=l ^ ' 



Vr G 



Now, if we take a time-reparametrisation t = t{T) such that 



it follows that Xt\ , Xt^ G V and, consequently, Xrl — Xt-^ G V. Taking into account 
that 

1=1 ^ ' 
we have X*]"^ - Xt^ = ^"^^ v'd/dv' = G 1^. 
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On the other hand, consider a family of parametrizations t = tx^ti{T~), with A G 
and ti G M, satisfying the conditions 



Consequently, the family of vector fields 

i=i ^ ^ 

is included in V. Note that, for every ti, the above family of vector fields can be considered 
as a curve in V. As is a vector space, all the derivatives of such a curve, that is, the 
vector fields 

HP dP f) 

— \X^^ (x, v)] = y —{X^F\t,, X, \-'v))— = XVx, v), 
d\p ^ nv n d\p^ ^ " dv' P^ ' 

i=l 

are included in V . □ 



6 Superposition rules and systems of HODEs 

In order to introduce a general theory of superposition rules for systems of HODEs, let 
us recall some basic concepts of the theory of higher-order tangent bundles p3T]. 

Given two curves p, a : R — !■ M" such that p(0) = cr(0) = Xq G M", we say that they 
have a contact of order s at xq, with s G N, if they satisfy 

^^(0) = ^^(0), j = l,...,s, 

for every function / G C°°(M"). The relation 'to have a contact of order s at Xq' is an 
equivalence relation. Each equivalence class, say t^^, is called an s-tangent vector at xq. 
Now, we define T^^R" as the set of all s-tangent vectors at xq and we put 



U 



- xq'- 



It can be proved that (T^R", vr, M"), with vr : t^^ G T^M'" ^-^ Xq G M", can be endowed 
with a differential structure of fibre bundle. Let us briefly analyse this fact. 

Every global coordinate system {x^, . . . , x^} on R" induces a natural coordinate system 
on the space T^'R". Indeed, consider again a curve p. The s-tangent vector t^^^ associated 
with this curve admits a representative 

t do^ t^ d^o^ 
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which can be characterized by its coefficients 

In consequence, the mapping : t^^ G T'^R" h-)- (xq, ?/o^^*, • • • , ?/o'^'**) G R('^+^)" gives a 
canonical global coordinate for T'^M". Obviously, the map vr becomes a smooth submersion 
which makes T^'M" into a fibre bundle with base R". We hereby denote each element of 
T^R" by = 

Now, given a curve c : t G R i— ?■ c{t) G R", we call prolongation to T'^R" of c the curve 
: t G R t*c(t) G T'^R'^, associating with every to the corresponding equivalence class 
of c(t — to), and given in coordinates by 

Definition 6.1. We say that a system of s-order ordinary differential equations on R" 
given by 

^ = F^(tx- ^] 1 = 1 n (61) 

dtp V '^^^'■■■'t^^^-v' ' 

admits a superposition rule if there exists a map T : (T*~-'^(R"))''" x R**" — )■ R" of the form 

X = T(t^^^j, . . . , ^xi^^y ^1; • • • ; ^sn), 

such that the general solution x{t) of system (16. ip can be written as 

x{t) = T (t^-ia;(i)(t), . . . , t^-'x(„)(t); h, . . . , Kn) , (6.2) 

for any generic family of particular solutions X(i)(t), . . . ,X(rn)(t) of (16. ip and /ci, . . . , ksn 
being a set of constants related to the initial conditions of each particular solution. 

Note 6.2. Observe that, according to the above definitions, we have T^R" ~ TR" and 
t^a;(a)(t) = [x(^a)(t),dx(a)it)/dty Hence, Definition 13.31 describing superposition rules for 
systems of SODEs turns out to be a particular case of the above definition. Moreover, if 
we put T^R" = R", Definition 16 . 1 1 reduces to the standard superposition rule notion. 

Definition 6.3. We say that a system of ordinary differential equations (16. ip of the order 
s is a HODE Lie system if the first-order system 

dt ^ ' 

dt ^ ' i = l,...,n, (6.3) 

^^ = F^(t,x,y«,...,y(^-)), 

obtained from (16. ip by adding the new variables y*^-'^* = d^x'^/dP, with i = 1, . . . ,n and 
j = 1, . . . , s — 1, is a Lie system. 
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Observe that the results and definitions described in previous sections can be directly 
generalized to systems of HODEs. This is why, instead of detailing such generalisations, 
we shall merely describe a simple but relevant result ensuring the existence of superpo- 
sition rules for HODE Lie systems. In next sections, this result is used to determine a 
superposition rule for second- and third-order Kummer-Schwarz equations. 

Proposition 6.4. Every HODE Lie system admits a superposition rule. 

Proof. Note that every solution of system (16.31) is of the form t^~^Xp{t) for a particular 
solution Xp(t) of system (16.11) and vice versa. Consequently, the superposition rule $ : 
(T'-^M")™ X W -> T"-^M" for (Q allows us to write that the general solution t'-^x(t) 
of (16. 3 P in the form 

t'-^x{t) = $(t^-^a;(i)(t), . . . , t^-ix(„)(t); h, . . . , 

in terms of generic families of particular solutions X(i)(t), . . . , X(m)(t) of (16. ip . their deriva- 
tives up to s — 1 order, and the constants ki,...,ksn. Applying the projection vr : 
rps-ij^n _^ iQQf^]^ sides of the above relation, it follows that the general solution 

x{t) of system (16. ip can be written as 

x{t) = (tt o $)(t"-ix(i)(t), . . . , t^-^X(^)(t); fci, . . . , ksn). 

In other words, T = tt o $ : ('j'^^-'^M'*)™' x — )■ M*^ is a superposition rule for (16. ip . □ 



7 Examples of superposition rules for systems of SODEs 

Let us illustrate now the results described in the previous sections by means of various 
examples extracted from the physics and mathematics literature. As the first simple 
instance, consider the n-dimensional isotropic harmonic oscillator 

^ = -u\t)x\ z = l,...,n, (7.1) 

with a time-dependent frequency u{t). This system appears broadly in the physics litera- 
ture. For instance, it occurs in the study of the fluctuations of the tachyon field obtained 
by using effective Lagrangians [321 [33], in the description of the movement of a particle 
on a heated spring |35, in the analysis of the properties of diverse interesting nonlinear 
differential equations, like Milne-Pinney equations, with many applications in physics 

[HIiniEs]. 

As shown in Section [HI systems (17. ip admit a base-superposition rule. Consequently, 
Proposition 15.41 ensures that the systems of the form (17. ip must be SODE Lie systems. 
Actually, this can be proved easily. In view of Definition 15. H demonstrating that each 
system of the form (17. ip is a SODE Lie system reduces to proving that every first-order 
system 

dx'^ 

^ = l,...,n, (7.2) 

dV n, , 
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is a Lie system. Any system of the above type describes integral curves of the time- 
dependent vector field Xt = Xi + uj'^{t)X^, with 

i=l i=l ^ ^ i=\ 

spanning a Lie algebra Vho of vector fields isomorphic to s[(2, M) [15]. It turns out that 
every system fl7.2p is a Lie system and, therefore, the isotropic harmonic oscillators with 
time-dependent frequency (17. ip are SODE Lie systems, as Proposition 17. 101 states. 

As equations (17. ip admit the same base-superposition rule. Proposition 15.61 ensures 
the existence of a finite-dimensional Lie algebra including Vuo and = 
Indeed, it is a straightforward computation to check that A^, X\, X2, and X3 generate a 
Lie algebra of vector fields isomorphic to Qi{2, M). 

Let us now turn to exemplifying that superposition rules for systems of SODEs need 
not be invariant under time-reparametrisations, as pointed out in Section |H Recall that 
Milne-Pinney equations (I3.10p admit a quasi-base superposition rule depending on a 
constant of motion I^. A time-reparametrisation t = Jq g^^^ ^dt' transforms (I3.10p into 
the dissipative Milne-Pinney equation 



^2 jrp J -2F{t{T)) 

±1 = _^(t(^))e-^(*W)^ _ u;2(t(^))e-2F(t(r))^ ^ ^ (7 3) 

and the superposition rule (13. lip for Milne-Pinney equations yields that the general 
solution of the above equation can be cast in the form (13. lip again, but in terms of a new 
constant 1 3 reading 



e 



This proves that the (quasi-base) superposition rule for Milne-Pinney equations is not 
invariant under time-reparametrisations. Moreover, as SODEs (17.30 admit, generically, 
a time-dependent superposition rule (see [ID] for details), they need not be SODE Lie 
systems. Indeed, it was proved in [2] that systems (17. 3p are not SODE Lie systems. 
Then, from Lemma 15.51 Milne-Pinney equations cannot admit a base-superposition rule. 

Let us now derive a new superposition rule for a relevant type of (nonautonomous) 
second-order differential equation: the second-order Kummer-Schwarz equation 

^ = 7^f5V-2&oX=^ + 2ao(t)a:, (7.4) 



dt"^ 2x \ dt J 

where 60 is a constant and ao(t) is an arbitrary time-dependent function. The study of 
these, hereafter KS-2 equations, is motivated by its appearance in the theory of super- 
position rules [T7] and in the analysis of second-order differential equations, where they 
appear related to the so-called Kummer problem [IQ]. These equations also appear asso- 
ciated with the so-called second-order Gambler equation [22] and can be used to describe 
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certain cosmological problems [H] . Moreover, the solution of several cases of KS-2 equa- 
tions amounts us to solving certain Milne-Pinney and Riccati equations [2n |22| [i2]. As 
these equations are ubiquitous in the physical literature, e.g. they appear in cosmology, 
quantum mechanics, classical mechanics [HlllS], the study of KS-2 equations can be con- 
sider as a useful approach to the analysis of these equations and their respective related 
physical problems. 

In order to describe a superposition rule for KS-2 equations, we shall first prove that 
these equations are SODE Lie systems, which ensures, by Proposition [5]2l that they admit 
a superposition rule and indicates how to derive it. 

Recall that demonstrating that KS-2 equations are SODE Lie systems relies on proving 
that the first-order system 



dx 






= V, 


~dt 


dv 


_ 


H 


~ 2 X 



(7.5) 

2box^ + 2ao{t)x, 
is a Lie system. To do this, consider the vector fields 

Xi = 2x^, X, = x^ + 2v-^, X, = v^+(^--2box')-^. (7.6) 
ov ox ov ox \2 X J ov 

Since 

[Xi,X2] = Xi, [Xi,X3] = 2X2, [X2,X3] = X3, 

they span a Lie algebra of vector fields isomorphic to s[(2,M) and, as system (17. 5p is 
determined by the time-dependent vector field 

^* " ^ i 27 ~ ^^""^^ ^ 2aoit)x\ — = X3 + ao(t)Xi, 

the second-order Kummer-Schwarz equations are SODE Lie systems. 

It is interesting that, like time-dependent frequency harmonic oscillators, Kummer- 
Schwarz equations are SODE Lie systems related to a Vessiot-Guldberg Lie algebra iso- 
morphic to s[(2,M). This can be used to establish interesting relations between these 
equations and other (SODE) Lie systems associated with the same Lie algebra (cf. [T^). 

Once it has been proved that KS-2 equations are SODE Lie systems, the following 
step toward deriving their superposition rule is, in view of Lemma 15. 2[ to determine the 
part of the standard superposition rule for system (17.51) describing the x coordinate of its 
general solution. To do this, let us apply the method described in Section [2j 

The vector fields Xi,X2,X3 form a basis for a Vessiot-Guldberg Lie algebra of (17. 5p 
and their prolongations to (TM)^ are linearly independent at a generic point. Let Xi, X2, 

and X3 be diagonal prolongations to (TM)^. As [Xi,X3] = [Xi,X3] = 2X2, if a function 
F : (TM)^ -> M satisfies XiF = X3F = 0, then X2F = 0. Thus, obtaining a common 
first-integral for Xi,X2,X3 reduces to finding a common first-integral for Xi and X3. 
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Consider the canonical coordinates {a;o, fo, f i, X2, ^2} in (TM)^ and suppose that 
the common first-integral F for Xi and X3 depends only on the variables xq^xi^vq, and 
Vi. As XiF = 0, the method of characteristics yields that F must be constant along the 
solutions of the characteristic system 



Xo 



dvi 

Xi 



dxo = dxi = 0. 



(7.7) 



Integrating the above system, we find that that there exists a certain function F2 : — ?■ 
M such that F (xq, xi, fo, fi) = F2 {xo,xi,^ = xiVq — xqVi). In terms of the variables 
xq, Xi, ^, f 1, the condition X^F = X3F2 = reads 



^ + ViXo\ dF2 , dF2 , 
+ ■^1^^ \- 



Xi 



dx(] 



dxi 



3 fe^r2^vix^\ 

XiXq ) 



+ 2ho{x\xQ - xlxi] 



0. 



As F2 does not depend on vi, the above equation implies 



^dF2 

Xi dxo 



+ 



2xiXo 



+ 2bo{xlxo - xlxi] 



dFj, 



xodF2 ^ dF2 ^ 3^dF2 
xi dxo dxi xi dC, 



Applying again the method of characteristics to the second equation, we see that there 
exists a function F3 : — )■ M such that F2{xo,xi,$,) = F^i^Ki = xo/xi,K2 = Let 
us express the first of the above equations using the variables Ki,K2, and ^. As a result, 
it turns out that 



x\ 



dKi 



2Ki 



+ 2boKl{K^ - Kf] 



dFs 



0. 



The characteristic system for the preceding equation is 

dK2 



dKi 



^^+2b,Kl{K,-Kl) 



(7.8) 



whose solution is 



{vqXi 



ViXq) 



3 3 

1 



^2 _|_ ^2 
XoXi 



Hence, Fi is a first-integral common to X\ and X3. Similarly, if we suppose that F depends 
only on xo,a;2,fo, and f2, or, alternatively, on a;i,X2,fi, and f2, two new first-integrals 
appear: 



{vqX2 -V2XQf xl+xl 

- + 40o . 

X0X2 



XqX2 



T,= ^^^^^l^^ + 4b,^-i±^. (7.9) 

" " X1X2 



/y>3 /y»3 

XIX2 



Since d{Ti,T2)/d{xo,Vo) 7^ at a generic point of (TM)^, the procedure described in 
Section |2] allows us to determine the values of Xq and Vq in terms of Xi,X2,fi,f2, and 
two constants giving rise to a superposition rule for system (17. 5p . Indeed, fixing Fi = ki 
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enables us to determine the value of vq in terms of xo,xi,vi, and ki. Substituting this 
value into the equation k2 = T2 and with the aid of T^, we can express xq in terms of 
xi,X2, ki, k2, and Fs as 

^ _ {T^ki - 8bok2)xi +{rsk2 - 8boki)x2 ± 2Afc„fc, (13)^3X1X2 - 4bo{xl + xl)]'/^ 
^° 166or3 + x^^X2'[{kiXi - k2X2y - 64bl{xl + xl)] 

(7.10) 

where 

Xk,Mi^3) = [25Qbl + k,k2T, - Aboikl + kl + Vl)]^'^ . 

Expressions (17. 9p ensure that [r3XiX2 — 46o(a;f + 2:2)]^'^^ is real. Meanwhile, for each pair 
of solutions Xi()f:), X2(t), ki and ^2 niust be chosen so that XkiM^^^) ^^^1- 

Since r3 depends on the variables Xi,X2,fi,f2, it is clear that expression (I7.10p con- 
stitutes a part of a superposition rule for any system of the form (17.51) . describing the 
component x of its general solution in terms of two particular solutions Xi(t), X2{t) of 
(17.41) . their derivatives fi(t),U2(t), and two constants ki and /c2. Therefore, in view of 
Lemma this allows us to write the general solution x{t) of the equation (17.41) in terms 
of two particular solutions, their derivatives, and two constants. This provides us with 
a superposition rule T : (t^.^, t^^; fci, ^^2) G (TR)^ x x = T(ti.^, t^^; fci, A;2)) G R for 

KS-2 equations of the form 

^ (Jfci - 8&ofc2)xi + {Ik2 - 86oA:i)x2 - 2Afc„fc,(/)[/xiX2 - 4&o(xf + xfiY''^ 
l^b^I + Xi^X2^[{kiXi — k2X2y — 646g(xf + x|)] 

where / = r3 is regarded as a function of the variables of (TR)^. Note in addition that 
the above expression can also be naturally considered as a quasi-base superposition rule 
of the form G(xi, X2, /; /ci, /C2) for KS-2 equations. 



8 A superposition rule for the third-order Kummer— 
Schwarz equations 

The present section is devoted to the study of third-order Kummer-Schwarz equations 
[171 Sa is] of the form 

with 60 being a constant and ao(t) being any time-dependent function. Our aim is to 
exemplify how the results of Section [6] can be applied to investigate a relevant third- 
order differential equation. As a result, it is shown that third-order Kummer-Schwarz 
equations, hereby KS-3, are HODE Lie systems, and an interesting relation to KS-2, 
Riccati, and Milne-Pinney equations is pointed out. Finally, a new superposition rule for 
KS-3 equations depending on a single particular solution and three constants is derived. 
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The relevance of the study of KS-3 equations rehes, for instance, on their relation to 
the so-called Kummer's problem [T71 HH HS], Milne-Pinney equations [16], and Riccati 
equations [HI [21], |16] . These relations can be used to study multiple physical systems 
described by these latter equations through KS-3 equations, e.g. the case of quantum 
non-equilibrium dynamics of many-body systems [H]. Furthermore, Kummer-Schwarz 
equations with 60 = can be rewritten as {x,t} = 2ao(t), where {x,t} is the so-called 
Schwarzian derivative of the function x{t) with respect to t [48]. 

In order to study KS-3 equations, let us define y^^^ = dx/dt, y^'^^ = d'^x/dt'^, and write 
equation (18.11) in the form 



dx 
11 




dy^'^ 
dt 


= y''\ 


dy(^^ 


3y(2)2 


, dt 


2yW 



(8.2) 



2W^' + 2ao(t)yW. 



(1) 



Consider the vector fields Xi, X2, and X3 on T^ 



X2 = y 



dy(^y " dyW ^ dy(^y 

d 



(1)3' 



d 



.3) 



satisfying the commutation relations 

[Xi,X3] =2X2, [X2,X3] 



x. 



[Xi,X2] =Xi. 



Obviously, these vector fields span a Lie algebra of vector fields isomorphic to s[(2,M). 
Additionally, system (18. 2p describes integral curves of the time-dependent vector field 
Xt = X3 + ao(t)Xi. Thus, KS-3 equations are HODE Lie systems. Let us derive a 
superposition rule for them. 

The vector fields Xi, X2, X3 are linearly independent at a generic point of T^M. There- 
fore, the diagonal prolongations Xi,X2,X3 of Xi,X2,X3 to (T^M)^ span a generalized 
distribution P. Such a generalized distribution is three-dimensional in a neighbourhood 
of a generic point, where the distribution becomes regular. Hence, the vector fields of the 
distribution admit, at least locally, three common first-integrals. As [Xi,X3] = 2X2, we 
have [Xi,X3] = 2X2, and obtaining first-integrals common for all the vector fields of V 
reduces to determining first-integrals common for Xi and X3. 

Let us first analyse first-integrals of the vector field Xi on (T^M)^, i.e. solutions 



F : (T2 



of the equation 



(1) 



dF 



(1) 



dF 



dy\ 



(2) 



0. 
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The method of characteristics shows that the first-integrals of the above vector field are 
functions constant along the solutions of the so-called characteristic system of Xi, namely, 

— = — T^r, dy^^^ = dy[^^ = dxo = dxi = 0. 
Vo Vi 

Such solutions are given, in an implicit form, by the algebraic equations ^ = Vo^yi^ — 
yi^yo \ Vo = yo \ vi = yi \ yo = xo, yi = xi, where ^,vo,vi,yo,yi are certain real 
constants. In other words, any first-integral F : (T^M)^ — )■ R of the vector field Xi 
depends only on the previous variables. Hence, there exists a function F2 : ^ M such 
that F{xQ,xi,y^Q\yf^ ,y^Q\yf^) = F2{yo,yi,vo,vi,^). 

Remember that we are interested in determining a common first- integral for the vector 
fields Xi and X^. In view of the above result, = amounts us to 

,(2) 



where we defined ai = y\ ' . In terms of the vector fields 

Zi = + vi- — + h 2bo{Vf)Vi - v^vo) ^, 

dyo dyi vx Ovq \ 2viVo ) di 

OVq OVi ot, 

we can write X3F2 = Z1F2 + ^^2-^2 = 0. Since F2 does not depend on ai, the previous 
decomposition implies Z1F2 = Z2F2 = 0. Applying the method of characteristics to 
Z2F2 = 0, we find that F2 must be constant along solutions of the characteristic system 

dvo dvi 1 1 n 

= = ayo = dyi = 0. 

Vo vi 3^ 

Therefore, F2 depends only on the variables Ki — vi/vq, K2 = vf/^, yo,yi, i.e. there 
exists a function F3 : — )■ R such that F{xo, Xi, vq, Vi, ao,ai) = F^lyo, yi, Ki, K2). 

To obtain a common first-integral for all the vector fields in D, it remains to impose 
ZiF — Z1F3 = 0. In the coordinate system {yo, yi, Ki, K2, ^, ai}, this equation reads 



.i/3^v3 (±m m Kim 

^ ' \K,dyo dyi K2 dK, 



= 0. 



The characteristic system corresponding to the above equation is 

%i K2dKi KfdK2 



dyo 



Ki Kf \Kf-2hoKl{l-Kl)' 
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From the last equality of the above system, we obtain 
dK2 



dKi 2Ki Kl 



Kl 



V^iri-46o(i + i^f; 



for a certain real constant Fi. Consequently, a common first-integral of the vector fields 
of V reads 

Kl + Ah^Kl{l + Kl) 



Now, dy\ = K2dKi/ Kf and the above expression yields 



lA) 



where sg stands for the sign function. Assume, for simplicity, sg(i^2) = 1 and bo < 0. 
Hence, 

T2 = (^-86oA'i + Fi + A^Abl{l + Kf)-boK,T,^ e-2xiv^ 

is a second first-integral. Likewise, from Fi and expression dyo = K2dKi/Kf, one gets 

dyo _ 
dKi ~ 



K.JT^K,- Aboil + Kl 



I.e. 



yo 



2V^o 



In 



-860 + i^iFi + A^Abl{l + Kf)-boK^T^ 



is another first- integral. As 9(Fi, F2, T3)/d{xo,yl)^\ yj^^) 7^ at a generic point of (T^M)^, 
fixing /ci = Fi, and k2 = F2, we can easily express Ki in terms of ki, ^2 and Xi. Using 
this and putting ^3 = F3, we obtain 



xo = k3 + In 



2V^om - fl,^{x,)] 



6Abl{ki - 2e2v^^'iA;2) - kifl^f^^{xi) + 8bo{6Abl - kf + e^v^^ifci 



-bo 



where fk^Mi^^) = ki — e^^^^~^^^k2. Note that the above expression is a superposition rule 
for third-order Kummer-Schwarz equations, which provides the general solution xo{t) of 
any instance of such equations in terms of a generic particular solution Xi{t) and the 
constants ki,k2,k3. Obviously, this represents an improvement with respect to other 
similar expressions for KS-3 equations, which allows us to describe their general solutions 
in terms of two particular solutions of a time-dependent frequency harmonic oscillator 
[T7] . In addition, this expression is an instance of a quasi-base superposition rule for a 
third-order differential equation. 
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9 Conclusions and Outlook 



We have proposed and analysed a general concept of a superposition rule for systems 
of HODEs. Some specific types of such superposition rules that appear in the litera- 
ture have been studied and other new types have been introduced and investigated. All 
our results have been illustrated with examples extracted from the mathematical and 
physics literature. In particular, two new superposition rules for second- and third-order 
Kummer-Schwarz equations have been derived. 

There are still many open questions concerning the properties of superposition rules for 
systems of HODEs. For instance, it would be interesting to find methods for analysing 
the existence of solutions of system (14.11) . which would facilitate the determination of 
the existence of superposition rules for systems of SODEs. Additionally, it would be 
interesting to apply the methods developed here to analyse first-order systems from a 
new perspective. 

In the future, we intend to study the whole Riccati hierarchy [49J, some of whose 
members, like second-order Riccati equations, have already been analysed by means of 
the theory of Lie systems p!0]. Further, we aim to apply our results in the analysis 
of soliton solutions of PDEs described by the Riccati hierarchy [H]. Additionally, we 
plan to employ the theory of Lie systems so as to geometrically explain the relation 
of Kummer-Schwarz equations to Lie systems associated with a Vessiot-Guldberg Lie 
algebra isomorphic to s[(2,M). This may be used to clarify their known connections with 
time-dependent harmonic oscillators or Riccati equations [2n HI] as well as to establish 
new ones. These and other topics will be analysed in forthcoming works. 
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